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Galerkin  Method  for  Solving  of  Singular  Integral 
Equation  of  Diffraction  Problem* 

Smirnov  Yu.G.,  Tsupak  A. A. 


1 The  statement  of  the  diffraction  problem 


Let  P = {i : 0 < £i  < a,  0 < X2  < b,  0 < x3  < c}  be  a resonator  with  perfectly  conducting  boundary.  Let 
Q be  a three-dimensional  body,  located  in  P.  Q is  characterized  by  tensor  permittivity  e and  constant 
permeability  /z 0.  We  suppose  that  components  of  £ are  smooth  functions  in  Q and  is  invertible 

in  Q;  QDdP  — 0.  Let  P/Q  be  homogeneous  and  isotropic  medium.  Incident  and  diffraction  fields  depend 
on  time  variable  as  e~tut. 

We  will  find  electromagnetic  diffraction  fields  E and  H,  satisfying  Maxwell’s  equations  in  P \ dQ  : 

rot  H = - iojiE  + jfj.  , .. . 

rot  E = iuj/xH  - jft  . 


The  complete  field  should  have  continuous  tangent  components  at  dQ  : 


and  must  satisfy  the  following  boundary  condition: 

E'\qp  = 0 . 


(2) 


2 Integro-differential  equations  for  the  diffraction  problem 

We  will  express  the  solution  of  the  stated  problem  in  terms  of  vector  potentials  AE  and  Ah  [4]: 

Ae  = J GE{x,y)JE(y)dy  ,AH  = J GH{x,y)]H{y)dy  , 

Q Q 

E — iu>y,0AE  - j^grad  div  AE  - rot  Ah  , 

H = iojeoAn  - grad  div  Ah  + rot  AE  ■ 


(3) 


Here  Je=J%  + Je,  Jh  =Jh+3h’  He*  3h  are  polarization  currents).  GE,GH  are  Green  functions  for 
Helmholtz  equation,  conforming  to  the  arbitrary  currents  JeJh- 

Ge,Gh  are  known  [3]  to  have  the  form  of  diagonal  tensors  (the  components  of  GE  are  written  out 
below): 

ftl  V 2 f„  . .firn  \„:nZ7rm  \ /irn  t f slvyX3  sll7(c  — y3),  X3  < y3 

°E  C0i(  « Xl)  sm(  6 X2) C0b^  o 2/1)  bin(  6 2/2)|sh7y3  sh7(c  - x3),  X3  > y3 


n=l  m=0 


G%=  E £ 


n=  1 m=l 


air)' shy  c 


sin(^xi)  sin(^a;2)  sin(^f  t/i)  sin  (™y2) 


sh7X3  sh7(c-y3),x3  < y3 
sh7y3  sh7(c-x3),  x3  > y3 

ch7x3  ch7(c  - y3),  x3  < y3 
ch7y3  ch7(c  - x3),  x3  > y3 


(4) 
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Here  7 = \J (~-)2  + (- j 1)2  — (the  proper  branch  for  square  root  is  chosen  as  in  [2],  §2.3),  £q  = 1 and 
en  = 2 for  n = 1,2,3,.. . . 

We  can  obtain  the  following  integro-differential  equations  (under  the  condition  p = /j0I  in  P ): 


£(x)  = E°(x)  + k20  J Ge  ^ - 1 E{y)dy  + grad  div  j GE 


e(y) 
£ 0 


-I 


3{y)dy , 


and  we  have 


H(x)  = H°(x)  — iu>£0rot  [ Ge  — / E(y)dy,  x € Q 
J Eo  J 


(5) 


We  can  extract  singularity  of  Green  function  G.  Using  Fourier  transformation  and  interpolation  poly- 
nomials we  can  obtain: 

l etko|3— »l  „ 

Ge(%, y)  = r ■ I + diag{gx{x,y),g2{x,y),g3(x,y)}, 


Ar  \x  — y | 


where  gu  are  smooth  functions. 


3 Galerkin  method 

Let  us  introduce  the  following  auxilary  function 


- OO  OO 

G(*,V)  = -E  £ ^sin(^x1)sin(^x2)sin(^2/1)sin(^2/2)x 

n=l  77i=l  /g\ 

x f shqx3  shq(c  - y3),  x3  < y3  ; 

\sh7y3  sh7(c  - x3),  x3  > y3 


(7) 


The  derivatives  of  G are  connected  to  the  derivatives  of  G'E  through  the  equalities: 

dGl  dG  . , „ „ 

dXi  ~ dyS 

Before  describing  the  method  itself  we  should  make  some  transformations  of  equation  (5).  Denoting 
- ij  as  £ and  (sM  - fj  £ as  .7  we  obtain  the  following  equation 

AJ  :=  |j(x)  -kl  J GEJ{y)dy  - grad  div  J GEJ{y)dy  = E0(x)  (8) 

Q Q 

We  can  write  vector  equation  (8)  as  a system  of  three  scalar  equations: 

^^liPix)  - kl  j GlE{x,y)Jl(y)dy  - ~divx  j G{x,y)J{y)dy  = El0{x),  1 = 1, 2, 3.  (9) 

i=1  Q ‘ Q 

We  will  determine  the  components  of  approximate  solution  J in  the  following  way: 


N 


N 


N 


j1  = Y,akfk(x),  j2  = £)  &*/*(*)»  j3  = 5]cfc/f(x), 


(10) 


k= 1 


fc=i 


k=  1 


where  flk  are  basis  ” hat” -functions  dependent  essentially  on  x'.  The  explicit  form  of  f\  is  given  below. 

Let  Q be  a parallelepiped:  Q = {x  : ai  < x1  < a2,bi  < x2  < b2,ci  < x3  < c2},  Q c P.  We  will  cover 
Q with  smaller  parallelepipeds 


njfcim  = = 4-i  < xl  < 4+1-  4 < z2  < 4+1'  xm  < x3  < x3m+1} 

1 <Z2  Gti  o , _ 62  ~ bl  . >>  Co  ~ Cl 

= ax  + — - — fc,  x\  = bi  + 2 — - — l,  xzm  = ci  + 2—^ — m; 


ai) 
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where  k — 1, . . . ,n  — 1;  l,m  = 0, 1, . . . , | — 1. 

Denoting  ( x*,  — %k- 1)  as  hl  we  get  the  formulas  for  /£Im  : 


/ klm  ^ 


-xfe-i 


, if  x1  € [x^;  x£]  and  x € 11^ 


klm 


*T+i- 


or 


r^r>  i/  x1  € [x^ ; x£+1]  and  x € n^(m 
l 0,  i/x^ni>im 


_ J i/xen[Im 

1 0,  if  x $ njIm 


Functions  f£lrn  and  f%lm  should  be  determined  by  similar  formulas.  Since 

/klmlx’ €{xj[_,  ,xP,}  — /fc(mlx2£{xf_,,xj  } — 0>  /fclmlxs6{i®i_1,x^ 


,}=0’ 


(12) 


(13) 

(14) 


every  component  of  approximate  vector  solution  vanishes  at  some  side  of  Q.  However  the  constructed  set 
of  basis  functions  does  satisfy  the  necessary  approximation  condition. 

Introducing  total  enumeration  for  basis  functions  we  get 


flflfb  k = l,...,N, 

where  N — \{n3  - n2). 

It  is  convenient  to  represent  the  augmented  matrix  for  determining  unknown  coefficients  ak,bk,Ck  in 
block  form: 


( 

A12 

A13 

Bl  \ 

I A2i 

A22 

a23 

Bx 

(15) 

\ a31 

A3  2 

A33 

B:  ) 

where  columns  Bk  and  matrices  A*;  are  determined  by  formulas: 


B'k  = {Eq,  fi);  (16) 


A:  — 1,2,3;  i — l,. . . ,N.  (f,g)  determines  the  scalar  product  in  L2,  (/, g)  — f f(x)g{x)dx. 

Q 

Applying  the  formulas  of  integration  by  parts  to  both  internal  and  external  integrals  and  taking  into 
account  (7)  and  (14)  we  obtain: 


Akt=  j £kiflj{x)fi(x)dx  — Skik2  J J G%(x, y)fj {y)/* (x)dydx— 

n(nnt  nt  n'. 

-j S °^y)^[y)^{x)dydx- 

rP  n‘. 

* J 


(18) 
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